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!*> ■ Abstract 

O ,' It is demonstrated, how the theoretical analysis preceding the nu- 

merical calculations helps to calculate the energy of the ground state of 
^NJ \ helium atom, and enables to avoid qualitative errors in calculations of 

the characteristics of double photoionization. 

o 

O ■ 1 Introduction 

As it stands now, many of the publications devoted to interaction of photons 
and electrons with atoms and molecules have similar structure. In the first step 
a most general quantum mechanics formula for the cross section is written. In 
the next step the most "accurate" numerical functions for the bound state 
electrons are found. Since the binding energy can be determined with high 
accuracy, the wave functions pass the test for reproducing the binding energies 
with a good accuracy. In the next step computer is at work. Sometimes there 
is a possibility to try the influence of the final state interactions. The usual 
choices are: plane waves, Coulomb functions... 

About twenty years ago such approach was justified to some extent. Most 
calculations were connected with the characteristics, which could be detected 
in the experiments. In the real experiments the targets obtained the energies 
of the order of the electron binding energies. There was no small parameter, 
and thus no possibility to evaluate the equations, determined by the original 



formalism of quantum mechanics. The part of atomic physics, related to the 
real experiments was indeed becoming the science about the precise compu- 
tations. The near-threshold behavior of the processes, where one could find a 
small parameter, was rather an exception. 

The situation changed at the end of 80-th, when new synchrotron sources 
of the photons became available. The experiments with the photon energies of 
about 10 keV went on. For the light atoms these energies are much larger then 
the binding energies. A small parameter emerged in the theoretical problems, 
which could be of practical importance. 

The general approach to the "high energy atomic physics" can be based 
on the works of Bethe jT]. A somewhat modernized presentation is given in 
j2 E] • Being initially formulated for the electron-atomic scattering, these prin- 
ciples can be expanded to any process of the interaction between the high 
energy particles with the bound systems. The wording "high energy" means 
that the energy exceeds strongly the binding energy of the system. The cross- 
sections of such processes can be expressed through certain parameters of the 
bound systems. The interactions between the fast and slow participants can 
be treated as perturbation. Each act of transferring of large momentum leads 
to a parametrically small factor. In the region, where the processes are kine- 
matically allowed for the free electrons, a small momentum q of the order of 
the average momentum of the bound electron 77 is transferred to the nucleus. 
The processes outside this kinematical region ( which we shall refer to as Bethe 
surface) require large momentum q ^> rj to be transferred to the nucleus. The 
processes outside the Bethe surface are thus strongly quenched. 

Hence, the amplitudes of the high energy atomic physics processes depend 
on the parameter 

k = rj/q . (1) 

On the Bethe surface k is of the order of unity, while k<C 1 outside the Bethe 
surface. Using the system of units with h — c — 1, we can write for the lowest 
state of the hydrogenlike atom 

rj = maZ , (2) 

with m being the electron mass, a = 1/137 is the fine structure constant, while 
Z is the charge of the nucleus. We shall use this relation for estimations in the 
many-electron cases as well. 

The perturbative treatment of the fast electrons, which is another point of 
the Bethe theory, leads to certain misunderstandings sometimes. The sections 



"Perturbation theory" of the textbooks on quantum mechanics often treat 
only the short-ranged forces. Indeed, iteration series term of the Lippman- 
Schwinger equation for the electron wave function in the Coulomb field of the 
nucleus contains the infrared-divergent terms. In this sense, there is no such 
thing as "perturbative series for the Coulomb field" and the words " the fast 
electrons can be described by the plane waves " become confusing. However, 
the solution was found long ago. It was guessed by Dalitz jl] and proved by 
Gorshkov [S], that the divergent terms compose the factor exp(i&) with real <J> 
in the amplitudes, thus cancelling in the cross-sections. Hence a way to solve 
the problem is to present the electromagnetic field of the nucleus in the form 
V(r) = —aZe~ Xr jr with A — » 0. Now the Lippman- Schwinger equation and 
the perturbative theory become legitimate, while the terms which depend on 
A (i.e. InX) cancel in the cross-sections. 

The perturbation parameter of the nonrelativistic electron with the three- 
dimensional momentum p, moving in the Coulomb field is jH] 

iz = maZ/p . (3) 

The condition £# <C 1 means that the relative contribution provided by this 
interaction to the electron wave function is much smaller than unity jHj. Since 
different regions of distances may be important for different terms of the ex- 
pansion of the wave function, one should be careful and expand the amplitude 
in powers of £z- The small distances are connected with transferring of large 
momentum q ^> rj to the nucleus. This process can be considered as pertur- 
bation as well, with k 4 <^ 1 being the parameter of expansion-Eq(l). 

However, one often finds large deviations between the calculations with 
the plane waves and Coulomb field functions applied for the description of the 
continuum electrons, even at large values of energies. Sometimes this leads 
to the statement that "plane wave approximation never works". However, 
as explained in |7j, there are two parameters, which depend on the electron 
momentum p. Besides the parameter £§ which is the true parameter of the 
expansion of the wave function at the distances of the order of the size of 
the atom, there is a parameter 7r£z> which enters through the normalization 
factor of the Coulomb wave function anyway. In the processes with the large 
momentum transferred to the nucleus, there are some more terms, containing 
the parameter 7r£z- Fortunately, they also compose a factor. Thus, assuming 

< Z « i, e z « i (4) 

one can carry out expansion in powers of £f . Also, the factors depending on 
the parameter n£z cancel in the ratios of the cross sections of the processes 



containing the high energy electrons with the close values of energies. The 
recent calculations, presented in [S] illustrated the general validity of the latter 
statement. 

The attempts to include the final state interactions of the outgoing electrons 
in a straightforward way also met some difficulties, caused by the long-range 
nature of the interaction. The approach which enabled to include the e — e 
interaction in the lowest order of the parameter of their interaction 

e = {maf/pl (5) 

with p r being the linear momentum of the relative motion was worked out in 

Since the high-energy experiments were not available in atomic physics 
for a long time, these ideas were not applied much. The only group which 
systematically used and developed the approach was that of V. G. Gorshkov 
and his collaborators jTUj, who focused on the Coulomb field calculations. 

Now one can conclude that it is not enough to require that the bound state 
wave function provides an accurate energy value. It should provide the ade- 
quate description of the electron density on the distances, which are important 
for the considered process. If the process is not allowed for the free electrons, 
the important distances can be of minor importance for determination of the 
binding energy. In treating the final states one can use the perturbative expan- 
sions in powers of £§, £ 2 and k 2 , taking care that all the terms of the leading 
order are picked in the amplitude. 

It looks to me, that the story of investigation of the double photoionization 
of helium provides good illustrations of how the application of this approach 
enabled to obtain the results. It illustrates also, how ignoring of these ideas 
led to the erroneous results. The theory gave two lessons: at the end of 60-th 
and at the end of 90-th. The relativistic aspects of the problems will not be 
touched. I will show , how the main results can be obtained by the standard 
methods of nonrelativistic quantum mechanics. I will show also, that the 
Lippman- Schwinger equation appears to be a powerful tool. However, I start 
with the example of a static problem of the calculation of the ground state 
binding energy of the helium atom. 



2 Calculation of the ground state energy of 
helium atom 

The binding energy of the ground state of helium atom is known from the 
experiments with very high precision. The relative accuracy was 10~ 6 in late 
50-th [2. It is about 2 ■ 1CT 7 today [TT]. This stimulated the attempts to solve 
the Schrodinger equation for the helium ground state as precisely as possible. 

The Hartree-Fock (HF) method, developed in early 30-th, reproduces the 
value of the binding energy with the relative error of about 1.5%. The experi- 
ments are much more precise. The reasons of this failure are known. The HF 
function does not depend on the distance between the electrons. The proba- 
bility to find the two electrons to be close to each other is thus overestimated 
in HF approach. 

An alternative approach was developed in the same years as the HF method. 
It consists in employing approximate wave functions ip a { r i-, r 2-, r i2)) presented 
by explicit analytical formulas. Here ri,r 2 denote the distances between the 
electron and the nucleus, r 12 is the distance between the electrons. The lower 
index a stands for "approximate". These functions contain a number of pa- 
rameters. The values of the parameters are chosen to minimize the expectation 
value 

E = (iP a \H\iP a ) , (6) 

of the Hamiltonian H of the helium atom. The simplest example is the product 
of two hydrogenlike functions in the field of the nucleus with certain "effec- 
tive nuclear charge" Z e . The effective charge plays the role of the variation 
parameter. In this case 

i>a(ri,r 2 ) =4>h(r i)ip h (r 2 ), (7) 

with the well known analytic expression for iph{r) = e~ ar , a = maZ e . ( All 
the wave functions are normalized to ip(0) = l).The solution Z e = 27/16 [S] 
provides the energy value with the error of 2%. 

More complicated functions of the form 

Mn,r 2 ) = e- a ^ +r ^P( ri ,r 2 ,r 12 ) (8) 

with a and the coefficients of the polynomial 

P(r h r 2 , r 12 ) = Y,c ijk r\r 3 2 r\ 2 (9) 



being the variation parameters, have been used since late 20-th. The earlier 
versions contained 3 parameters. By the end of 50-th the functions with more 
than 200 parameters were employed [2|. The functions P, containing also the 
dependence on the ratios , r ] 2 , and ( ri ~p> have been used. This enabled to 
obtain the convergence of the results to a certain value, interpreted as the true 
energy value, with the precision of several units of 10~ 10 , reproducing 5 signs 
in the value of the binding energy ( The calculations with the accuracy better 
then 10 -4 require the account of the finite mass and motion of the nucleus, 
inclusion of relativistic corrections, etc). 

However, Barlett et al. [12] showed that the functions (8), with P having a 
polynomial structure could not satisfy the Schrodinger equation for the ground 
state of helium. Such solution does not exist since it does not have a proper 
limit at the triple coalescence point r% = r 2 = r\i = 0. The functions which 
treated the three-particle singularity of the Hamiltonian were found by Fock 
|13j . They took the form 

P(n, r 2 , r 12 ) = T,c ijkn r[rir k 12 [ln{rl + r 2 ) 2 ]", (10) 

thus containing the logarithmic terms. (Usually, the Fock function is presented 
by using the hyperspherical variables |14j). 

Inclusion of the logarithmic terms has little effect on the energy value. 
However, account of these terms improves the convergence of the variation 
calculations. Also, inclusion of such terms enables to diminish the number 
of variation parameters. For example, the calculation with the accuracy 10~ 9 
required 1078 parameters of the function (9), while it was sufficient to use the 
function (10) with 52 parameters to achieve the same precision [T^j. Further- 
more, as noted by Myers et al. [HI], for some of the functions (9) the variation 
procedure chooses the parameters in such a way, that the expansion mimicked 
the function (10) with a smaller number of parameters. 

Behavior of the true ground state wave function \I/(ri,r2, r^) at the two- 
particle coalescence points is determined by the Kato cusp conditions ^7j. In 
the point of electron-nucleus coalescence r± — it is 

r 5 = - Z n°> r 2 5 r 2), (11) 

or 1 

while in the electron-electron coalescence point rn = it is 

^ = ^,01, (12) 



with r = — standing for the Bohr radius. Equations (11,12) can be viewed 
as the conditions of the cancellation of the singular terms in the Schrodinger 
equation. 

Let us illustrate this, considering the single-electron wave function ip(r) of 
the bound electron in Is state of certain effective field U(r), which approx- 
imates the interactions with the nucleus and with the other electrons. The 
Schrodinger equation for this case is 



U)"(r) + -V>'(r)) + U(r)ip(r) = Eip(r). (13) 



-1 /,«, , 2 

1m 



At r — > the interaction U(r) is determined by interaction with the nucleus. 
Thus, U(r — > 0) ~ — — . Hence, for r->0we can write 

^' {T)+ ^(±m-!*M) =E Hr). d4) 

Since both sides of Eq. (14) should be regular, the expression in brackets 
should have a finite value. Thus 

ro Mr) = - zm , (i5) 

at r — > 0. Similar analysis of the Schrodinger equation for the atom of helium 
leads to Eqs.(ll),(12). Note, that these relations could have been obtained 
immediately after the Schrodinger equation was written in 1926. Surprisingly, 
they were found only thirty years later. 

The proper treatment of the two-particle singularity is even more important 
than the account of the three-particle one (THl- This is because the three 
particle coalescence point is felt in the smaller volume of integration of the 
matrix element (6). Myers et al. [TH] suggested the simplest form of the wave 
function, which satisfies both Kato conditions 

Mn,r 2 ,r 12 ) = e -a(n+r 2 )+6m (16) 

with a = maZ, ft = 22. It is amusing that this function, which does not con- 
tain fitting parameters , provides the same precision in determination of the 
binding energy, as the function (7), which contains a variation parameter Z e . 
The function (16) knows nothing about the three-particle singularity. The sim- 
plest function, which would account both two- and three-particle singularities 
is described by a more complicated formula, also presented in [IB| . 



If the approximate wave function is chosen in a proper form, the variation 
procedure exhibits the tendency to find the parameters, providing small errors 
in the Kato conditions (11,12). The wave functions 

Mri, T2, r l2 ) = e- a ^ + ^ +br ^c ljk r\rir k 12 , (17) 

(i + j + k < N) analyzed by Teng and Shakeshaft [IH] provide the relative 
errors of about 6% and 25% in Eqs.(ll, 12) correspondingly, even in the case 
of 4 parameters (N=l). Although the Kato conditions are satisfied only ap- 
proximately, it helps in improving the convergence. The function (17) with 14 
parameters (N=4) provides the value of the binding energy with the relative 
deviation 3 • 10 -5 from the one, obtained with the functions, described by Eqs 
(8,9) and containing 210 parameters [2|. 

In the nowadays studies the subject of investigation is usually not only the 
expectation value (6), but also the local energy functional 

E(ri,r 2 ) = 



<i'a 



which gives more detailed information about the relation of the approximate 
function ip a to the true wave function \l/. Another improvement consists in 
developing the approaches, based on direct solution of the Schrodinger equa- 
tion. In one of the methods, suggested by Haftel and Mandelzweig ( see [T^j , 
where the other methods are also discussed) the wave function is presented as 
a product of a correlation factor, describing the coalescence points singulari- 
ties and a smooth function, expanded in a series of hyperspherical harmonic 
functions. The detailed description of these points is beyond the scope of this 
lecture. 

I hope that this Section made one convinced that taking care of the proper 
analytical structure of the approximate wave functions, one obtains the better 
results. Now we shall see that this becomes increasingly important for the 
dynamical problems. 



3 Asymptotics of the double photoionization 
cross section: a; -7 / 2 or uj~ 5 ^ 2 ? 

The double photoionization of the ground state of helium is studied since the 
late 50-th. The process provides a direct probe of three-body physics, since the 

8 



three charged particles are involved. A huge number of experimental results 
gives a good possibility for theoreticians to test their approaches. 

The mechanism, which determines the main contribution to the nonrela- 
tivistic asymptotics of the cross section a 2+ (u) with u standing for the photon 
energy, was clarified by Kabir and Salpiter [20]. At the values of u, exceeding 
strongly the binding energy /, the process can be viewed as a single photoion- 
ization, followed by the transition of the second electron to continuum. The 
latter process takes place due to the change of the Hamiltonian of the system. 
This is the analog of the shake-off (vstryahivanie) , suggested by Feinberg [21] 
for the description of ionization during the nuclear beta decay, due to the sud- 
den change of the electromagnetic field, caused by the change of the charge 
of the nucleus. However, using the same wording shake-off, one must remem- 
ber, that, in contrast to the beta decay, the mechanism can not be treated as 
due to the sudden change of certain effective field. Byron and Joachain |22] 
showed that, being attributed to the sudden change of self-consistent (Hartree- 
Fock) field, the mechanism underestimates the asymptotic value by a factor of 
3. The initial state correlations beyond the effective field are thus of crucial 
importance. 

The electron, which interacts with the photon directly, absorbs most of the 
energy t\ ~ u). The secondary electron carries the energy of the order of the 
binding energy €2 ~ I <C u>. Thus, the energy dependence of the amplitude 
and of the cross section of the high energy nonrelativistic photoionization are 
the same, as that for the single photoionization process. 

Recall now the mechanism of single ionization. 

3.1 Asymptotics of the single photoionization 

Here we assume that the electrons are described by the single-particle wave 
functions. The general expression for the amplitude of ionization of Is state is 

F + = (4vra) 1/2 / d 3 riP*(r)-fe t{kr) iP(r), (19) 



with ip p i being the continuum wave function with the asymptotic momentum 
Pi, if) describes the bound Is electron, while 

7=-<^, (20) 

m 



is the operator of interaction between the photon and electron, e stands for 
the photon polarization. The Fourier transform provides 

F+ = ^af 2 J^r Fpl (q)^M<l-k), (21) 

Here the functions ipp are the Fourier transforms of the functions ip. The 
process requires that momentum 

Q = k-p x (22) 

is transferred to the recoil ion. Note, that the energy of the fast electron 
e i = Pi/2m ^> /, while the binding energy I — %—, with 77 standing for the 
average momentum of the bound electron. Thus, e±^$> I means, that p\ ~^> 77, 
i.e. the momentum of the fast electron is much larger than that of the bound 
electron. Also, the energy dependence of the photon and electron momenta are 
k = uj and p\ ~ (2mw) 1 ' 2 . Hence, Q 3> 77 in this process. The large momentum 
(i.e. momentum, exceeding strongly the average value for the bound electron) 
should be transferred to the recoil ion. 

Start with the plane wave 

TpF P i(q) = ^f ] p M = (27r) 3 5(p 1 - q), (23) 

for the description of the final state electron. In this case all momentum Q is 
transferred to the recoil ion by the initial state electron. The contribution to 
the amplitude is 

F+ = (AnaY^-^MPi) ■ (24) 

m 

For the hydrogenlike functions the asymptotics ?phF(Pi) = — — t^ h 0(pf 6 ) 

at pi ^> 77 can be obtained explicitly. One can expect the p^ 4 behavior to 
be true in the general case, since the momenta p\ ^> i] correspond to small 
distances r <^C 77 -1 , where an electron function ip has a hydrogenlike shape. 
Anyway, this can be proved rigorously, and the pre-asymptotic coefficient can 
be calculated. Since the s state wave function does not depend on the direction 
of the vector r, we can evaluate 

^j F (p) = />r?Kr)e-^ = ~% [ drr^(r)\cos(pr)]' = ... - ^itM + 0{v~ % ) 
J p 2 J p 4 

(25) 
Here the last equality is obtained after using two integrations by parts. 
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To get convinced, that Eq.(25) indeed presents the asymptotics of the am- 
plitude, let us try to include the interactions between the outgoing electron and 
the recoil ion. The lowest order perturbative correction to the outgoing elec- 
tron wave function, caused by interaction with the nucleus V^can be obtained 
by using the standard formula of quantum mechanics 

v c ?M°) 

with the sum over the states m, which do not coincides with n, e„ )m stand for 
the unperturbed energy values. This provides 

with the contribution 

, . ,-,/, f d 3 q (eq) SiraZm , , . , , 

^ = -^ ' I Wf-^r w> - a»(ri - ?) Mq) ' (28) 

to the amplitude. In the amplitude F± both initial and final state can transfer 
large momenta q or p\ — q to the nucleus. The small factor of the order p^ 4 
is partially compensated by the large phase volume of integration of the order 
p\ , multiplied by the factor maZ. Thus, the amplitude F^ is indeed £z times 
smaller then F + , dropping faster with the increasing photon energy. 

Thus, corrections to asymptotics provide the terms 0(£z), dropping as 
ijj 1 ! 2 . The large numerical coefficient n can be understood, e.g. by considering 
the procedure of normalization of the function i^ p i(r) = i/jL '(r) + ip^ (r). 

Calculating the value V P i(0) = ^?(0) + Vff(O) = / -§^[^f pX {q) + Vgi(g)] 
one obtains: 

V P i(0) = l+i z J (29) 

with 

oo 

J = -/-ln^- = vr/2, (30) 

' ' 

and thus ip p i = 1 + 7r£z/2. Additional power of n can be viewed as coming 
from —iln(-l) = it if the integral J is evaluated in the complex plane. Hence, 
the amplitude _F + provides the asymptotics of the total amplitude. This leads 
to the well known asymptotics of the cross section [2S] 

a + (u) = duo- 712 (31) 
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Because of the large corrections of the order ir£z, the nonrelativistic asymp- 
totics can be achieved only for very light atoms. In other cases the relativis- 
tic treatment becomes important at the energies far below the nonrelativistic 
asymptotic region. 

3.2 Asymptotics of double photoionization. Spurious 
controversy 

The general expression for the double photoionization amplitude can be written 



as 



F 2+ = (*/| 7 |*i>, (32) 



with \l/j and *£?f standing for the exact wave functions of the initial and final 
states, while 

7 = 7i + 72; lj = * — (33) 

m 

describes the interactions of the photon with the two-electron system. The 

asymptotics of the amplitude (2QJ is determined by the interactions in the 

initial state, which transfer small energy e 2 ~ I to the secondary electron. 

Large momentum 

Q = k-p,-p 2 (34) 

is transferred to the nucleus. 

To obtain the asymptotics, we can describe the final state by the product 
of the single-particle functions, putting ^/j(r 1 ,r 2 ) = Vpi v 1)^/^2(^2 )• Here ijj c 
stands for the continuum wave function in the Coulomb field of the nucleus, 
while ip^ denotes its lowest order of expansion in powers of £z (the plane 
wave). Since p\ 3> P2, we can put 7 = 71. 

Analysis, similar to the one, carried out in the previous subsection, provides 
F 2 + = (Ana)^-^ f *L*1%®*fI&, & (35) 

The integral on the rhs is saturated by small q ~ r\ <C p\- The calculation, 
similar to the one presented by (25), provides ^p(pi, q) ~ Pi 4 at Pi 3> V- Thus, 
at P2 ~ f] the dependence of F 2+ on the photon energy is the same as that of 
F + . The integral over the energy distribution is saturated by e 2 ~ I and thus 
does not depend on the photon energy. Hence, the double photoionization 
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cross section a 2+ obtains the same energy dependence as the single ionization 

a 2+ (uj) = C 2 uj- 7/2 (36) 



cross section, i.e. 



In paper J22] the double photoionization cross section was calculated by 
employing the product of the Coulomb functions for the description of the 
final state. The asymptotics was obtained by expansion of the amplitude in 
powers of £#. It was noted that the result depends on the form, assumed for 
the operator of interaction between the photon and electrons. For the exact 
solutions of the Schrodinger equation the velocity form, presented by Eq.(33) 
is equivalent to the length form 

7 = 7i + 725 lj = Wfii&Tj) (37) 

with ujfi being the difference of the energies of the final and initial states, w/j = 
uj in our case. While the velocity form provided the asymptotics presented by 
Eq.(36), the length form gave 

a 2+ (uj) = C£uj- 5/2 . (38) 

To understand this, note that since p 2 < Pi, the important space region is 
r i "C r 2 . Thus we can put 7 = iuj(er 2 ). Unlike the single photoionization, 
the electron, which interacts with the photon directly, does not obtain angular 
momentum. The double photoionization amplitude takes the form 



(4™)^ / _|_^^ (gi) ^ ;(g2)(eg2)v]/F(gi5g2) _ (3Q) 



d 3 q 1 d 3 q 2 

If we describe the fast outgoing electron by the plane wave - Eq.(23), the 
dependence of the amplitude on the momentum of the fast electron pi ~ 
(2771a;) 1 ' 2 is determined by the factor 

LU^ F (p 1 ,q 2 ) ~ ^ > ( 40 ) 

Pi 

containing additional energy dependent factor u; 1//2 , compared to Eq.(24). This 
leads to Eq.(38) for the asymptotics of the cross section. In the calculations 
of [22] such behavior remained true, when the fast electron was described by 
the Coulomb field function. 

The authors of 22 j suggested that the to~ 5 ^ 2 result is spurious, and that 
the corresponding contribution should vanish if the exact wave function are 
used. This guess was proved by Aberg several years later |24j . 
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3.3 Solution of the controversy 

Here I present a simple derivation of the main results of J21] . Large momentum 
Q « pi can be transferred to the nucleus either by the initial state electron, 
or by the final state electron. In the length form each of these mechanisms, 
being treated separately, provides the contribution of the same order in the 
parametric expansion. The wave function of the outgoing electron can be 
described by the superposition of the plane wave and the lowest order Coulomb 
correction. Electron interactions in the final state as well as the higher order 
Coulomb corrections contribute beyond the asymptotics. The two terms of the 
expansion of the final state wave function provide the two contributions to the 
amplitude F 2 L + = F|+ + F|+, i.e. 



d 3 q 2 

Of 



Fit = (4™) 1/2 / ^r Fp M)m)A k {q 2 ) W 



with 



M&) = "j ^^i*(?i)^(gi, Q 2 ) (42) 

and k = 0,1. For the plane wave contribution we find immediately 

A (q 2 ) = uj^ F (p u q 2 ) , (43) 

with the large momentum p\ 3> rj being transferred to the nucleus by the 
initial state electrons. (Recall that the lower index F denotes the Fourier 
transform). The amplitude A is determined by large qx — Pi ^> r\. The 
integral, corresponding to A\ is saturated by small momenta q± ~ rj. This can 
be seen immediately from Eq.(27), which presents the function i/j f L explicitly. 
The momentum q± can be neglected everywhere, except the wave function typ- 
This provides 

A 1 (q 2 ) = -u—± \ — ^ F {qi, g 2 ) = -^pf(0, q 2 ). (44) 

Pi J [Zir) p 1 

Here we introduced the partial Fourier transform 

^PF(r,q 2 ) = /d 3 r^(r,f 2 )e-^ 2) . (45) 

The calculation, similar to that, presented by Eq.(25), gives 



stvoj 



A (q 2 ) = r * PF (0,<? 2 ). (46) 

p\ 
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Thus, the amplitudes F L q and F L f appeared to be of the same magnitude. 
(Recall, that this was not the case for the amplitudes Fq" and F x + due to the 
factor ^- on the rhs of Eq.(28)).This does not contradict to the fact, that 
the function ip^{r) is £z times smaller then ip^(r) in any space point. The 
amplitudes Ak(q 2 ), which can be presented as 

A k (q 2 ) = w/^ri^(fi)ttp F (fi,&), (47) 

are dominated by different regions of the space distances r±. While A is 
determined by small r\ ~ p^ 1 , much larger distances r\ ~ r/^ 1 provide the 
main contribution to A\. 

Treated separately, each of the terms Aqi would provide u;~ 5 / 2 law for the 
cross section. However, carrying out the inverse Fourier transform, one finds 
A + A\ = for any q 2 due to the Kato condition (11), and thus 

Fit + Fit = • (48) 

This shows, how the straightforward calculations in the length form, em- 
ploying the approximate wave functions, which do not satisfy Eq.(ll), provide 
the spurious tu~~ 5 ^ 2 asymptotics for the double photoionization cross section. 



3.4 The Lippman-Schwinger equation 

In this subsection I show that the Lippman-Schwinger equation (LSE) is a 
convenient tool for the asymptotical analysis. The general form of LSE for the 
two-electron system is 

^ = #(°) + G e V^ . (49) 

Here V stands for the interaction, e denotes the energy value, corresponding 
to the solution ^ of the Schrodinger equation, ty^ ' and G e are the solution 
and the Green function for V = (free motion). 

To simplify the notations, let us write the LSE for the single-particle wave 
function ip, which we assume to be the exact solution of the Schrodinger equa- 
tion in a certain effective field U(r). It is 

^ = ^(°) + G e Uip . (50) 

Now G e is the free single-particle Green function. In the momentum space 

(P\GV)= %~ f) M . "-0. (51) 

e — p z /2m + iv 
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For the bound states i[)^ = and 

^ = 2^/1^^®' < 52 » 

with e standing for the single-particle binding energy. In the asymptotics 
2me <C p 2 . The integral is saturated by q ~ rj. Thus, Eq.(52) presents the 
bound state wave function iPf(p) at large momenta p ^> rj in terms of the same 
wave function at small momenta q ~ 77. 

For s states we can put (p\U\q) = (p\U\0). Since large p correspond to small 
r, only the interactions with the nucleus are important in the asymptotics. 
Thus, putting (p\U\q) = — ^pr-, we obtain 

#p = — f^+oip 6 , 53 

Due to the Kato condition (15), this expression is identical to Eq (25). In 
similar way we can find the asymptotics of the wave function for the ground 
state of helium 

*f(P,? 2 ) = 2 + 0(p . 54 

pi 

with the partial Fourier transform \1/pf defined by Eq.(45). 

As to the wave function of the outgoing electron, the first iteration of the 
LSE provides the same result, as Eq.(26) of quantum mechanics (adding an 
imaginary infrared divergent term, which will be absorbed into the phase fac- 
tor, mentioned in the Introduction). Thus, calculating the double photoion- 
ization amplitude in the length form, we immediately find the cancellation, 
expressed by Eq.(48). Hence, in the LSE technique one would not obtain the 
spurious u^l 2 term in the asymptotics of the double photoionization cross 
section, even being not aware of the Kato conditions. 



4 Shape of the double photoionization 
spectrum: U or W ? 

4.1 Mechanisms of the process 

The progress of the experimental facilities in early 90-th made possible the 
measurements of the energy distributions of the double photoionization of he- 
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lium. The low energy measurements have been carried out. The results for the 
high energy are expected soon. This stimulated the theoretical investigations. 

Evaluation of the shape of the spectrum curve have been obtained in [25 j . 
The energy distribution da 2+ /de reaches its largest values at the edge of the 
spectrum, where £2 ~ I <C uj or e± ~ I <C uj. The shake-off mechanism is 
dominative here |20j. Large momentum Q ^> r\ is exchanged between the nu- 
cleus and the fast electron. Considering the larger values e\^ ~^> I, we need 
both electrons to carry large momenta Pi ^> 77. If the energies of the outgoing 
electrons are not too close, the process can be viewed as the single photoion- 
ization, followed by the electron scattering fZfi\ 127] . This final state scattering 
(FSS) mechanism dominates until we do not reach the very vicinity of the 
center of the spectrum, where the quasifree (QF) process becomes possible. 

To understand the point, recall that each act of the exchange by a large 
momentum leads to a small factor in the amplitude. The FSS mechanism 
requires the exchange by two large momenta. However in the vicinity of the 
center we can manage by a single large momentum exchange. The electron, 
which interacts with the photon directly, can exchange large momentum with 
the second electron without participation of the nucleus. Unlike the single 
photoionization, the process can go on with the free electrons. The Bethe 
surface is determined by the condition Q = 0, i.e. 

P1+P2 = k , (55) 

which requires 

5 = ^^ < x [^— « 1 . (56) 

uj y uj + m 

( Recall, that we do not consider very large values of the photon energies 
uj ~ m, which require the relativistic treatment of the outgoing electrons.) In 
other words, the free kinematics becomes possible in the vicinity of the center 
of the spectrum. 

Following the general ideas, presented in the Introduction, we expect the 
surplus of the spectrum curve near it's central point due to the described QF 
mechanism. The energy distribution is dominated by small values Q ~ rj of 
the momentum, transferred to the nucleus. 

To calculate the asymptotics, we can neglect the terms of the order Q/pi 
in the wave function of the outgoing electrons. The QF amplitude can be 
presented as 

Fqf(px,P2) = F QF (p,Q) = D(Q 2 )F fr (uj,5) (57) 
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with p = (pi —p 2 )/2. The factor D(Q 2 ) contains all the properties of the initial 
state. Its explicit form will be given below — Eq. (71). The amplitude Ff r of 
the absorption of the photon by two free electrons at rest can be presented as 

F fr (u,S) = (ep 1 )f(uj,S) + (ep 2 )f(u,-S), (58) 

since the electron system is a space-symmetric state. We do not need to clarify 
the form of the function f(u,S) here. Anyway, since (5< 1, the lowest order 
of expansion in powers of S gives 

F fr (u,0) = (e,n+p 2 )f(u,0) = (ek)f(uj,0) = 0, (59) 

with the lowest order nonvanishing term being 

F fr (u,S) = 2(ep)5f'(uj,0) . (60) 

Here /' denotes the derivative with respect to 5. 

The energy distribution in the vicinity of the center is the interplay of FSS 
and QF mechanisms. Their contributions change with the photon energy in 
different ways. At uj <C 2keV the FSS dominates, providing U shaped curve. 
At uj « 2keV the QF mechanism becomes more important. The curve obtains 
W shape. 

Note that usually the photoionization calculations are carried out in the 
lowest order of expansion in powers of kr c with r c being the characteristic 
value of the distances, which are important in the process. This is known as 
the dipole approximation. In the nonrelativistic photoionization r c ~ p^ 1 ~ 
(mu;) -1 / 2 and the dipole approximation corresponds to the lowest order of 
expansion in powers of k/p ~ (cu/m) 1 / 2 . Going beyond the dipole approxima- 
tion have been something unusual, at least until the recent time. The dipole 
approximation QF amplitude turns to zero due to Eq. (59). This happens 
because the two-electron system cannot carry the orbital momentum, equal 
to unity, in this kinematics. One should go beyond the dipole approximation 
to obtain the nonvanishing QF amplitude. Unfortunately, this point was not 
emphasized strongly enough neither in [25], nor in a more detailed paper [27j . 
The QF amplitude is expressed by Eqs. (57) and (60) with 

5 = M . ( 6 i) 

muj 



4.2 Dipole approximation calculations 

Meanwhile, several calculations of the high energy photoionization spectrum 
have been reported. They were based on the straightforward computation 
of the amplitude in the dipole approximation. The results appeared to be 
controversial. In [2H] the spectrum curve had U shaped form if the final state 
interactions (FSI) between the outgoing electrons were neglected. Account of 
the FSI provided W shaped curve. The authors of [22J faced the opposite 
situation. Their spectrum curve for to = 2.8keV had W shape with a very flat 
peak, when FSI have been neglected. The computations with FSI included 
provided U-shaped spectrum. Another series of calculations without FSI [3UJ 
EH] gave a peak in the center of the energy distribution. 



4.3 Explanation of the dipole approximation 
results 



To understand what happened J32J|33|, present the ground state wave function 
in variables R — [f\ + r 2 )/2 and p ' = f 1 — r 2 , introducing ty(R,p) = ^(fi,^)- 
Describing the outgoing electrons by the plane waves, we find for the amplitude 

F (0) = (4vr«) 1/2 I <?Rd i pe i( 8* t >-*V-'''fi)'yV(R, p). (62) 

Following the analysis, carried out in Sec. 3, we expect the asymptotics F^°'(p, Q) 
VG° 4 Q 4 ) at p, Q ^ r]. However, in the vicinity of the center of the energy dis- 
tribution we can make Q ~ rj. These values of Q will be dominative in the 
energy distribution. The calculation, similar to the one, presented by Eq. (25) 
provides 

F<°) = -(4™)^M . ^ f d s Reim) 99(R,P = 0) 

m p 4 J op 

at p ^> rj, Q ~ r), with d^(R, p = 0)/dp denoting the partial derivation of the 
function ty(R,p) with respect to p at p = 0. The energy dependence of the 
amplitude F^ is determined by the factor p -4 , caused by the large momentum 
exchange. 

Large momentum can be exchanged by the outgoing electrons as well. The 
lowest order perturbative correction ip^ to the final state wave function can 
be obtained by using the standard formula (26) with V c being replaced by 
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the interaction between the electrons. The corresponding amplitude F^ is 
of the same order of magnitude as F^°\ containing the same factor p~ 4 . The 
situation is similar to that with the shake-off amplitude in the length form, 
which was discussed above. The correction V G°) to the wave function of 
the outgoing electrons is £ times smaller than the free wave function ip^°'(p) 
in each space point p. However, the functions ip^ and ip^ contribute to the 
amplitude in different space regions. While the amplitude F^ is determined 
by small distances p ~ p~ l , the amplitude F^ is dominated by the values of p 
of the order of the size of the atom ~ 77 -1 . In both cases R ~ Q^ 1 ~ rj^ 1 , and 
both electrons are separated from the nucleus by the distances of the order of 
the size of the atom. 

Direct calculation provides 

F« = (47ra) 1/2 ^ • ^ / d 3 Re^ R ^(R, p = 0) (64) 

771 p 4 J 

at p ^> 7], Q ~ 7]. Each of the amplitudes F^ 0,1 \ treated separately, provide 
the contribution, which exceeds strongly the amplitude, describing the FSI 
mechanism and requiring Q ^$> 7]. Hence, each of the terms i^ 0,1 ), treated 
separately, would lead to a surplus in the energy distribution in the central 
region. However 



m v J 



m p 
x ^ W , = 0)-2n, a * (fi -; = 0) )=0, (65) 

since the last factor in the integrand turns to zero for all values of R due to 
Eq. (12). In the dipole approximation the leading contribution to the energy 
distribution is determined by Q 3> rj, providing the U-shaped spectrum curve. 
However, the true distribution can be obtained only by going beyond the dipole 
approximation. 

Now we can explain the results of |2H] — |HX] - The authors of J2H! used the 
ground state wave function with a weak dependence on p. It was too weak to 
provide a maximum, when neglecting FSI. Inclusion of FSI provided a spurious 
uncompensated W peak. In contrast, authors of [2H] used wave functions which 
satisfied the Kato condition (12). Thus they obtained U shape, when taking 
into account the FSI. This corresponds to Eq. (65). However, a spurious W 
peak emerged when the FSI were neglected. The spurious central surplus, 
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obtained in [3UJ EH] i n dipole approximation with FSI being neglected have the 
same origin. 

The Lippman-Schwinger equation (LSE), presented by Eq. (49) enables to 
obtain these results in a simple way. The amplitude F^ can be expressed in 
terms of asymptotics of the Fourier transformed wave function 

FW = -^a) 1 ' 2 ^ * F (-Q,p) , (66) 

m 

while the asymptotics of the function ^f(p, Q) is provided by LSE equation 

M-Q,fi = -—j^k *r{-Q,$ (67) 

at Q ~ r/, p ^> 77. The integral in the rhs is saturated by q ~ ?y. Hence, the LSE 
presented the function ^f(Q ~ r],p ^ rj) in terms of this function at Q ~ rj, 
q ~ rj. Using Eqs. (64) and (67) we immediately obtain the cancellation 
p(o) _|_ _p(i) — 0, without employing the Kato condition (12). 

Note that the interactions between the outgoing electrons appeared to be 
more important than their interactions with the nucleus. The latter provide 
only additional corrections, while the former determine a possible mechanism 
of the process. Similar analysis with the Coulomb functions for the continuum 
electrons as a starting point is given in J22] • 

4.4 Necessary properties of the approximate 
functions 

Most of the calculations of the energy spectrum are based on the original 
formalism of quantum mechanics. This means straightforward employing of 
the equation 

da = -\F\S(e 1 + e a - U -I + )t^7^ (68) 

with the amplitude F calculated as 

F = J d 3 ri d 3 r 2 r fa (n, f 2 ) ( 7 ie l ( fc>1 ) + 72 e i(fc>2) ) Mn, f 2 ) (69) 

by the direct calculation of the integral. One should inevitably use approximate 
wave functions ipf a and tf) a for the description of the final and initial states. 
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To reproduce essential physics, the approximate wave functions should posses 
certain properties of the true wave functions. 

As we have seen in Sec. 3, in the region of the shake-off domination e% m u, 
£2 ~ I, the FSI correlations are not important and one can approximate the 
final state by the product of the single-particle functions, putting ipf a (ri, f 2 ) = 
i J a Pl {ri)ip ap2 (r2) and employing the Coulomb field function for ip a p 2 - in the 
velocity form of the operator 7 one can use the plane wave for the description 
of the fast electron [20 . Description of the initial state by the product of the 
single-particle wave functions 

V>a(n,r 2 ) = i) s {ri)i) s (r 2 ) (70) 

leads to the quantitatively wrong results. Thus the initial state correlations 
beyond effective field approximation are important [22]. However, there is no 
formula, which would enable to select "good" wave functions. One should be 
even more careful, employing the length form of the operator 7. The lowest 
order Coulomb correction to the wave function of the fast outgoing electron 
should be included and Eq. (11) should be fulfilled for the initial state function 

Mn,r 2 ) [21. 

Turning to the part of the spectrum with e 12 3> I we see that the lowest 
order final state interactions should be included in the wave function ipf a . This 
is sufficient for the proper description of all the region of FSS domination. 
However, to describe the region of QF mechanism domination, one needs the 
wave function ip a , for which the Kato condition (12) is true. Otherwise the 
spectrum curve obtains a spurious W peak in the dipole approximation (the 
latter corresponds to replacement e lt - kri < 2 ^ by unity in Eq. (69)). 

As we saw earlier, the true W peak can be obtained beyond the dipole 
approximation. This is equivalent to evaluation of the rhs of Eq. (69), by 
including the second terms of the expansion of the exponential factors e l<ykrj ^ = 
1 + i(krj) + . . . . Such approach provides Eq. (57) for the QF amplitude with 
Ff r (u,5) defined by Eq. (60) and 

D(Q 2 ) = 2r [d 3 Re^ d ^ a{R ' P = 0) . (71) 

J dp 

The proper value of the derivative g' P ~ can be insured by the Kato con- 
dition (12). Thus, even taking care of separating the quadrupole terms, one 
needs the approximate function ip a , satisfying the Kato condition (12) for the 
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proper description of the QF central peak. Of course, any approximate func- 
tion with a nonzero value of ,q' P ~ provides a central peak describing the 
QF mechanism qualitatively. However, the quantitative results can be trusted 
only if Eq.(12) is satisfied. Considering such approximate function, which can 
be presented as combinations of the single-particle functions (70), we find the 
corresponding 

j> a (R, p) = 4> s (ti+ $ s (r - £j (72) 

to be the even function of p. Thus, the computations, carried out in the 
quadrupole approximation would show no trace of the QF contribution. The 
Hartree-Fock functions provide one of examples. 



5 Summary 

We saw that theoretical analysis, preceding the computer calculations, appears 
to be very helpful. In the static problem of the calculation of the ground 
state energy of helium atom it improves the convergence of the procedure and 
enables to diminish the number of phenomenological parameters. It helps to 
avoid "submarine rocks" in the calculations of the asymptotics of the double 
photoionization cross section. It enables also to find the true shape of the 
double photoionization spectrum curve. 

It was shown also that the Lippman-Schwinger equation is a powerful tool 
for investigation of the asymptotics. Unfortunately, the LSE and the Feynman 
diagram technique, which can be viewed as its visual illustration, did not yet 
become a part of everyday life it atomic theory studies. 

I am indebted to M.Ya. Amusia, A.I. Mikhailov, R.H. Pratt, T. Suric and 
especially to V.G. Gorshkov for helpful discussions. 
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